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1. Introduction Besides the conventional hadrons, quantum chromodynamics, the 
theory of strong interaction, predicts that the states containing gluonic excitations, 
namely glueballs and hybrid hadrons may exist. The confirmation of the existence 
of such states would be a strong support for QCD to be the true theory of strong 
interaction. Therefore, such states attract a lot of interest both from theories and 
experiments. Several experimental candidates of glueballs and hybrid states have been 
reported. For example, a candidate for the scalar glueball was reported by Crystal 
Barrel Collaboration at LEAR [1], and possible evidence for a 1 h light exotic hybrid 
candidate was reported by VES [2]. As experiments begin to focus on the search for 
glueballs and hybrid states, it's important to get more theoretical knowledge of such 
states from theory. In this paper we will concentrate on the hybrid states with heavy 
quarks. 

Many methods have been applied in the investigation of hybrid in the literature, 
which include the flux tube model [3], the bag model [4], the constituent gluon model 
[5], and lattice Monte Carlo simulations [6]. In bag model, the mass of the lowest 
hybrid mesons with light quarks is about 1.5GeV. The hybrid mesons with heavy 
quarks are studied by taking into account the bag deformation, and it's found that the 
masses of the lightest hybrid are 3.9GeV for ccg, and 10.5GeV for bbg. In the flux tube 
model, hybrid is interpreted as the phonon-type excitations of the string connecting 
the quark and antiquark. It predicts that the masses of the lowest hybrid mesons are 
1.8 - 2.0GeV for uug or ddg, 4.2 - A.5GeV for ccg, and \\.2GeV for big. 

The constituent gluon models for hybrid hadrons were introduced by Horn and 
Madula [5]. In these models the constituent gluon plays the same role just as the 
constituent quarks. It's assumed that the constituent gluon can have a finite effective 
mass m g and carries angular momentum l g . As a result, it predicts non-exotic quantum 
numbers for the ground states of hybrid. Exotic hybrid states including 1~ + and 0~~ 
are predicted in the higher- lying multiplet. 

In this paper, we study the hybrid mesons wiht heavy quark in the framework 
of B.S. equation. This method is also based on the idea of constituent gluon. We 
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only concentrate on the hybrid with heavy quarks. In contrast to the light quark 
sector, where the complicated flavor mixing between qq and qqg bound states may be 
important, and the qq spectrum itself is not very clear, for the heavy quark system, 
the QQ spectrum is well described by the the potential model, and the mixing between 
QQ and QQg is small. This implies that it is possible to find pure hybrid states in 
heavy quark sector. 

It's also interesting to note that in the constituent gluon model for hybrid, when 
the mass of the quark and antiquark tends to infinity, spin-symmetry would show up. 
We will make a more detailed discussion of this point later. 

The paper is arranged as follows. In section 2, we construct the B.S. equation for 
the hybrid states. In section 3, we construct the B.S. wave function of hybrid states. 
In section 4, the numerical result and a discussion are presented. 

2. The B.S. equation of hybrid states Let ip(xi) and ip(x 2 ) be the quark and an- 
tiquark fields at points x\ and X2, A^(x 3 ) the gluon field at point x 3 . Then the B.S. 
wave function is defined as 



X^i,x 2 ,x 3 ) = (0\T4>( Xl )4>(x 2 )A^x 3 )\M), (1) 

where \M) is the hybrid state with mass M and momentum P. The color indices has 
been suppressed. We define the relative and center of mass kinematic variables as 

x = + x 2 ) + V2X3, (2) 

x = Xi-x 2 , x = — x 3 , (3) 

2m Q m g 



where 771 = , rj 2 = 

2m Q + m g 2m Q + m g 

With the translation invariance, we have 

X^(xi,x 2 ,x 3 ) = e~ lBX Xix(x,x'). (4) 
We can further define the B.S. wave function in momentum space by the Fourier 
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transformation 

X p(P, q,k) = J d 4 xd 4 x'e^ x+kx '\,{x, x'), (5) 

where q is the relative momentum between the heavy quark and heavy antiquark, and 
k is the relative momentum between the constituent gluon and the center of the two 
heavy quarks. With a standard method we can obtain the B.S. equation for the hybrid 
state. 

1 if d 4 q' d 4 k' 1 

^ (P ' g ' k) = -h-m QV l-ml J J2nyj2^°^ ^ * > *' * > ^J^T Q 

(6) 

In obtaining the above equation, the propagators have been replaced by the free ones. 
Pi, P2 and ps are the momentum of quark, antiquark and gluon respectively, and we 
have the following relations 

p 2 = -P-q + ~, (7) 

P3 = V2P - k. 

The integral kernel in the above equation is defined as the sum of all three-particle 
irreducible graphs. However, it's impossible to calculate the kernel exactly from the 
first principle. In practice, some approximations and phenomenological assumptions 
have to be made in obtaining the kernel. As usual, we divide the kernel into two 
parts: the long distance part and the short distance Gfy. G® can be obtained 
approximately by calculating the following three one-gluon-exchange diagrams shown 
in fig. 1. 

Besides equation (7), we also have 

* = T p + ' + * 

p'^P-q'^ (8) 
Ps = V2P - k'. 
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Hybrid is color singlet. Each of the two-body system in the hybrid belongs to a definite 
representation of color SU(3). The Qg system belongs to a 3, the Qg system belongs to 

a 3*, and the QQ system belongs to a 8. The color factor for the Qg and Qg systems 

3 1 
is — -, while the color factor for QQ system is -. This implies that the force between 

the constituent gluon and quark or antiquark is attractive, while the force between the 

color-octet quark-ant iquark pair is repulsive and much week. Therefore, we have such 

a picture for the hybrid with heavy quarks: the constituent quark and antiquark are 

almost static with weak and repulsive force between them; the constituent gluon turns 

round the QQ pair and bind them together to form a bound state. Because the force 

between the quark and antiquark is very weak, in our treatment we simply neglect it. 

Then for the short distance part of the kernel, we have (see Fig. 1) 

GW(P, q, q>, k, k>) = (2tt) 4 5 4 ( Pi - ^^Lv^-f^i ~ ™ Q ) 

+ (2n)W(p 2 - p' 2 )^|^r^ 7 ^2 + m Q ), (9) 
where I — k — k', and T^ vp is the three-gluon-vertex, 

= (p 3 + P'z) P 9^ + Vz)v9ix P + P f 3 )^9u P - ( 10 ) 

As for the long distance part of the kernel, we have to construct it phenonemnologically. 
According to the experience with the ordinary mesons in the potential model, the long 
distance part of the kernel is mainly of scalar property. Therefore we only choose the 
first term of equation (10) as the spin dependence of the long distance part of the 
kernel, and we assume 

= -9pAP3 + Ps) ( (2vr) 4 5 4 ( Pl - p' 1 )(j/ 1 - m Q ) - (2vr) V(p 2 - p' 2 )(j/ 2 + m Q ) 

(11) 

where V = — , and — p— corresponds to a linear potential in the position space. 

Substituting the kernel into equation (6), and complete the integration over q', we 
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obtain 

f 3(4™,) + + p , } ^\ ^ ^ 



_L r d 4 k' ( 

<i-m Q J (2^ ^ 



(12) 

We see that now q is not a dynamical variable. By neglecting the q dependence of 
the propagators, this variable is integrated out. The propagator of the quark can be 
expressed as 

1 A+( Pl ) A-( Pl ) 

\Z ET + - — ^r)7o, (13) 



tfi-m Q pio - .Bi pio + #i 

1 d4) 



fa + rriQ p 20 + E 2 P20 ~ E 2 

where Ei = ^Jmf + p 2 , and A + (p) (A~(p)) is positive (negative) energy projector 
defined as 

A±(pi) = g.±7o(rP. + m Q ) ^ (18) 
For the constituent gluon propagator, we have 

pl-m] ~2E~ 3 { p 30 -E 3 ps + E 3 y (16) 
As an approximation, we neglect the negative energy part of the propagator. Such ap- 
proximation is reasonable especially for heavy particles. Then as usual, instantaneous 
approximation is made. Under this approximation we ignore the dk dependence of 
the kernel. After making k integration over the two sides of equation (12), we have 

2E 3 (M -E 1 -E 2 - E 3 )^(k) 

= " J j^y 3 [^y^T^Y - g^(Ps J 0^(A;')7oA-( P2 ) 

- A+( Pl ) 7 o J j^y 3 i^^T^Y + g^Ps + P^-V-^j (f) (k), (17) 
where 4>^{k) is the instantaneous three dimensional B.S. wave function. 

^(P,fc) = / dk oX ,(P,k) (18) 



3. The B.S. wave function of the hybrid In order to solve the B.S. equation, we 
should construct the B.S. wave function of the hybrid with the given quantum numbers. 
In the following we construct the wave function in heavy quark limit. The spin of the 
quark- ant iquark system, Sqq, can take two values, and 1. When Sqq = 0, we have 
the following hybrid states and its four dimensional wave functions. 
The 1+- hybrid: 

X ,(P, k) = (1 + Yh^e, + F 2 k-e k ±lt + F 3 k-eP,). (19) 

The 0— hybrid 

X ,(P, k) = (1 + fhsiF^ + F 2 P,). (20) 

The 1 hybrid 

X ,(P, k) = F(l + y^e^P'kie*. (21) 

The 2 hybrid 

Xli (P, k) = (1 + f^iF^k" + F 2VXl/ k x k»k ±p + F^k^P^). (22) 

When Sqq = 1 we have another series of hybrid states and its four dimensional 
wave functions. 
The ++ hybrid 

X ,(P, k) = (1 + Y)(F llfl + F 2 # ± k ±li + F4 ± P,). (23) 

The 1 ++ hybrid 

X„(P, k) = (1 + y){F 1 e livpal v e' , P a + Fis^pdWkl + Fss^evpok*^ 

+F 4 e^ pa YP p k CT k-e + F b e Xvpal x e l 'P»k° ' P p + F & e Xvpal x e» PWk^). (24) 
The 2++ hybrid 

X„(P, k) = (1 + Y)(FiV»vY + F 2 ri^k v # ± + Fw^VPr 
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+F A ri pa yk°k^ + F 5 rj pa kPk^ + F^W^^P^ + F^Wty^k^). (25) 
The -+ hybrid 

X M (P, k) = F(l + Y)e^l u P p k°. (26) 

The 1 h hybrid 

X M (P, fc) = (1 + y)(F lf /k ±lt + F 2 e-A;7 ±At + + F 4 ^ + F 5 $ ± k-ek ±li + F^k-ePj. 

(27) 

The 2~+ hybrid 

X M (P, fc) = (1 + f )(FiW^ a P^a + F 2 e flupa rj^P^ x + FsW^M^a- 

+F A e pup X X P p k a k x k ! 1 _ + F 5 e Xupa ri x P»yk° + F 6 e Xupa ri Xa P»yk°k a k ±p . 
+F 7 e Xupa r ] Xa P»yk°k a P p ). (28) 



In the above expressions e p is the polarization vector of the vector hybrid, and r] pv 
the polarization tensor of the tensor hybrid. q± p is a transverse vector defined as 
q± p = q p — q-VV p , and j± p is defined similarly, p's are scalar functions of q 2 , P-q. For 
simlicity, the same Fi in different hybrid states represent different functions. 

If we substitute the wave function into the B.S. equation, the spin symmetry can be 
seen explicitly. For example, for 1 ++ state, the components of its wave function, Fi, F 5 
and Fq are decoupled from the total wave function Xm an d the equations about these 
components are exactly the same as those of ++ hybrid. Similarly, for 2 ++ state, the 
components of the wave function, Fi, F 3 and F 4 are decoupled out from the total wave 
function, and the equations about these components are also the same as those of ++ 
hybrid. The reason is that when the interaction between QQ is week and negligible, 
the quark spin decouples from the light freedom in heavy quark limit. This result does 
not depend on the form of the kernel (10) and valid for more general case in which the 
interaction between QQ is omitted. Therefore, if the constituent gluon model is true, 
we wish to see this symmetry at least in the low-lying QQg hybrid. 
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4. The numerical results In the non-relativistic limit, in our numerical treatment 
we only choose the spatial part of the wave function. In the center of mass frame, the 
form of the wave functions are simplified significantly. 
For l + ~ state 

^(k) = (1 + 7o)7 5 (M + he-kh). (29) 



For state 



For 1 state 



For 2 state 



For ++ state 



For 1++ state 



(/>i(k) = fi(l + 7o)75^ 



<f>i{k) = A(l + lo)l5£ijkkje k . 



<t>i{k) = (1 + 1o)l5(fiVij k j + f2Vjkkjk k ki). 



(f)i(k) = (1 + 7o)(/i7* + hl-k h 



(30) 



(31) 



(32) 



(33) 



<f>i(k) = (1 + 7o)(/i£yfc7jefc + Sijkkjekk-^ + he ijk k^ k k-e + f^juk^ik^). (34) 
For 2 ++ state 



4>i(k) = (1 + 7o)(/i%7? + f2Vijkjk--f + f3Vjkkj-f k gi + UVjkkjk k -fi + f 5 rjjkkjk k k--fki 



(35) 



Substituting the above wave functions into the B.S. equation, we can calculate 
the spectrum of hybrid states. It's easy to see that the states l + ~ and (0, 1,2) ++ 
degenerate, and also the states (0, 1,2) and (0, 1, 2)~ + degenerate. 

It's easy to see that the long distance part of the kernel with the form of expression 
(11) is very singular at the point l 2 = 0. Some form of regularization is necessary. The 
method is to make the following replacement 

G2^-^(p3+P3)^f(27r)V(p 1 -p;)( ] / 1 -m Q )-(27r)V(p 2 -^)(^ 2 + m Q ) X " ' 



(l 2 + M 2 ) 2 




-((2 7 r)V(p 1 -p' 1 )(^ 1 -m Q ) 



(36) 



where u is a small quantity. In actual calculation, we let u — > 0. In this way, the 
infared divergence is subtracted out. 

In our calculation, the strong coupling constant a s is chosen as a running one. 



where the parameter a is introduced to avoid the infrared divergence. We choose 
a = 4.0 which implies that a s tends to its largest value 1 when I 2 — > 0. 

The other parameters including the string tension A, and the quark mass, m c and 
777.6, are determined by fitting the QQ spectrum, and we have: A = 0.2{GeV) 2 , m c = 
XAlGeV and m b = A.80GeV. 

The mass of the constituent gluon has been studied in a lattice calculation [7] , and 
in non-perturbative QCD based on gauge invariant Lagrangian [8]. These studies yield 
m g = (500 ~ 800) MeV. In our study, we adopt m g = QOOMeV. 

All the parameters are determined as above, we can solve the B.S. equation numer- 
ically. The results are list in the table I. 

As mentioned above, the choice m g = QOOMeV is reasonable. However, we can't 
say it's exact. We calculate the hybrid mass by changing m g in the the range of several 
hundred MeV . Fig. 2 shows the dependence of the mass of ccg on m g , and fig. 3 shows 
the dependence of the mass of bbg on m g . We can see that the hybrid mass is not 
too sensitive to the gluon mass. The mass of lightest ccg changes from 3.65GeV to 
3.80GeV when m g changes from 300MeV to 800MeV. In the same range of m g , the 
mass of the lightest bbg changes from 10.26GeV to 10.39GW. 

There is another kind of uncertainty due to the neglecting of the interaction be- 
tween the quark and anti-quark. Recall that for the ground state of charmonium or 
bottumonium in quark model, the interaction energy between the two quarks is less 



12tt 1 



27 ln(a + P/A 2 QCD ) 



(37) 
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than one hundred MeV, and here the strength of the interaction between the color- 
octet QQ is only one eighth of that between the singlet QQ. Therefore, we expect this 
uncertainty is small and can be neglected reasonably. 

5. Conclusion and discussion We have studied the hybrid states with heavy 
quarks in the framework of B.S. equation. In the heavy quark limit, we show clearly 
the existence of the spin symmetry which results in the degeneracy of the states of 
(0,1, 2)++ and also the degeneracy of (0, 1,2)~+ The states l+~ and (0,1, 2)++ arc 
also degenerate, and they are the ground state of hybrid, while the degenerate states 
(0, 1,2) and (0, 1,2)~ + belong to the first excited states. With a reasonable choice 
of m g (Q00MeV), we find that the ground state of ccg hybrid is near the threshold of 
DD meons, while the ground state of bbg hybrid is bellow the threshold of BB mesons. 
The energy gap between the first excited states and the ground states are 0.52GeV 
for ccg, and 0A8GeV for bbg, and such values are quite insensitive to the mass of the 
constituent gluon. 

In the experiment of e + e~ collision, several resonant states are found near 4.0GeV, 
such as -0(4040), 0(4160) and -0(4415). Within the conventional picture of charmonium, 
0(4040) is assigned as 3S state, and ■0(4415) is assigned as 4S* (or possibly 5S) state. 
But as has been argued by F. Close [9], 0(4160) perhaps is a 1 ccg hybrid state. 
In our model, the mass of the 1 ccg hybrid is 4.24GeV, so our result is roughly in 
agreement with this interpretation. 
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Figuer Caption 
Fig.l:The diagrams of B.S. kernel in the lowest order. 

Fig.2:The m g dependence of the mass of the ground and first excited ccg hybrid. 
Fig. 3: The m g dependence of the mass of the ground and first excited bbg hybrid. 
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Table Caption 
Table 1: The mass of the hybrid mesons (MeV). 
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Table 1. The mass of hybrid mesons(GeV). 



The ground state of ccg (1+- (0, 1, 2)++) 


3.72 


The first excited state of (ccg (0, 1, 2)—, (0, 1, 2)~+) 


4.24 


The ground state of (big l+~ (0, 1, 2)++) 


10.32 


The first excited state of (big (0, 1, 2) — , (0, 1, 2)-+) 


10.80 



